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ABSTRACT 

High resolution FUV echelle spectra showing absorption features arising from C I and CO gas in the 
(3 Pictoris circumstellar (CS) disk were obtained on 1997 December 6 and 19 using the Space Telescope 
Imaging Spectrograph (STIS). An unsaturated spin- forbidden line of C I at 1613.376 A not previously 
seen in spectra of /3 Pictoris was detected, allowing for an improved determination of the column density 
of C I at zero velocity relative to the star (the stable component), A( 3 P) = (2 — 4) x 10 16 cm~ 2 . Variable 
components with multiple velocities, which are the signatures of infalling bodies in the (3 Pictoris CS disk, 
are observed in the C I AA1561 and 1657 multiplets. Also seen for the first time were two lines arising 
from the metastable 1 D level of carbon, at 1931 A and 1463 A. The results of analysis of the CO A — X 
(0-0), (1-0), and (2-0) bands are presented, including the bands arising from 13 CO, with much better 
precision than has previously been possible due to the very high resolution provided by the STIS echelle 
gratings. Only stable CO gas is observed, with a column density iV(CO) = (6.3 ± 0.3) x 10 14 cm~ 2 . 
An unusual ratio of the column densities of 12 CO to 13 CO is found (R = 15 ± 2). The large difference 
between the column densities of C I and CO indicates that photodissociation of CO is not the primary 
source of C I gas in the disk, contrary to previous suggestion. 
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1. INTRODUCTION 



(3 Pictoris is the most extensively studied of the young 
planetary systems discovered in the last decade and a half. 
It is a bright Southern hemisphere star (type A5 V), lo- 
cated about 19.3 pc distant from the Sun, with a systemic 
radial velocity of 20 km s _1 (for a review of the (3 Pictoris 
system, see Vidal-Madjar, Lecavelier des Etangs, & Fer- 
let (1998)). It was observed in 1983 by the IRAS satellite 
to have a large excess of emission at infrared wavelengths. 
This was referred to as the Vega-like phenomenon and was 
identified as arising from an edge-on circumstellar (here- 
after CS) dust disk, presumed to be associated with plan- 
etary formation (Smith & Terrile 1984). It was soon de- 
termined through absorption spectroscopy that there was 
CS gas associated with the dust as well. 

A large body of evidence has accumulated indicating 
that there are comet-like bodies present in the (3 Pictoris 
CS disk. Collisions between dust particles are expected 
to produce submicron fragments which should be expelled 
from the system by radiation pressure on time scales much 
shorter than any plausible stellar age. Thus there must be 
a secondary source of particles; one model for the produc- 
tion of dust and gas in the CS disk focusses on evaporat- 



ing comets and is called the Orbiting Evaporating Bodies 
model (OEB). In this picture, the comets orbit the star 
at several tens of AU, and thus, the f3 Pictoris CS disk 
is a kind of "gigantic multi-cometary tail with its natu- 
ral constituents: gas and dust" (Lecavelier des Etangs, 
Vidal-Madjar, & Ferlet 1996). Spectra of (3 Pictoris show 
variable redshifted absorption features arising from gas in- 
falling toward the star at high velocities (and infrequently, 
blucshifted features as well); these features are best at- 
tributed to the evaporation of star-grazing comets, called 
the Falling Evaporating Bodies scenario (FEB) (Bcust et 
al. 1990). Also, gas at close to 20 km s _1 (zero radial ve- 
locity relative to the star) is identified in all observations 
and is called the stable gas component. This gas is difficult 
to understand, as modeling indicates that it should be ex- 
pelled from the system by radiation pressure; a continuous 
source for this gas is required. 

Neutral carbon and carbon monoxide have been ob- 
served in HST-GHRS UV absorption spectra of f3 Pic- 
toris; carbon monoxide is the only molecule detected in 
the CS disk to date (Vidal-Madjar et al. 1994; Jolly et al. 
1998). Since CO and C I can be dissociated and ionized 
by interstellar UV photons on time scales of the order of 
200 years, both must be continuously replenished. Carbon 
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monoxide, in particular, is difficult to reform after dissoci- 
ation in the f3 Pictoris environment. Thus, the presence of 
these species indicates that a secondary source for this gas 
should exist, just as for the CS dust (Vidal-Madjar et al. 
1994). Jolly et al. (1998) found the column densities of C I 
and CO to be comparable, around 10 15 cm~ 2 ; since their 
rates of destruction are also comparable, this was taken 
as evidence that the C I is produced by photodissociation 
of CO, which evaporates from comets orbiting at various 
distances and velocities. However, the C I column density 
was determined from a heavily saturated multiplct and 
is therefore quite uncertain. In the hopes of further con- 
straining the characteristics of C I and CO in the f3 Pictoris 
system, we have reinvestigated the transitions observed in 
the GHRS data, as well as some that were not seen due to 
the relatively low spectral resolution of GHRS compared 
with that of the STIS high resolution echelle. 

2. OBSERVATIONS 

HST STIS high resolution echelle spectra of (3 Pictoris 
were obtained on 1997 December 6 and 1997 December 19, 
covering the wavelength range from 1459 A to 2888 A in 
six exposures each day. Table 1 shows the log of observa- 
tions. All the absorption features discussed in this paper 
appear in either the first or second data set listed for each 
day (o4g001010/o4g002010or o4g001020/o4g002020).The 
data were initially reduced and calibrated using the STScI 
IRAF package calstis vl.8. Spectra with a signal-to- noise 
ratio of around 10 were achieved. Examination of the 
errors in the flux values showed that the error propaga- 
tion calculation had not been performed correctly and that 
thestated errors were too small. Therefore the data were 
re-calibrated using calstis vl.9a, correcting the underesti- 
mate of the measurement errors. 

The particular advantage of this data set over previ- 
ous comparable ones is the very high spectral resolution 
achieved. The instrumental line spread function using the 
E140H grating is well described by a Gaussian with a 
FWHM of 1 pixel, corresponding to FWHM = A/220,000 
(Sahu et al. 1999). However, since the detector under- 
samples the line spread function, the effective resolution 
to separate two adjacent lines with this grating is only 
R = 110,000. Using the E230H grating, the FWHM is 
about 2 pixels, corresponding to FWHM = A/110,000 
(Sahu et al. 1999). 

3. ANALYSIS 

3.1. CljXmSline 

This previously undetected spin-forbidden transition 
( 3 Po — ► 1 Pi) at 1613.376 A was clearly observed on both 
days. It is unsaturated, which allows for an improved de- 
termination of the column density of C I. As there was 
no change in the line between the two days of observa- 
tion, the spectra were averaged together to improve the 
S/N; the result is shown in Figure 1. Unfortunately, the 
other two fine structure lines in the multiplet were not 
reliably detected. Atomic data were taken from Morton 
(1991) and the continuum around the line was fit with 
a sixth-degree polynomial. Voigt line profiles were used 
to generate transmission functions, which were then con- 
volved with a Gaussian instrumental line spread function 
with FWHM = A/220,000 to create model spectra. The 



X 2 statistic between the model and the data was then min- 
imized to determine the best velocity centroid, v, column 
density of C I in the ground level, A^( 3 P ), and Dopplcr 
broadening parameter, b; error bars for these parameters 
were determined from the contours of x 2 . Since the line 
was unsaturated, the column density was also determined 
from the equivalent width of the line and b and v were de- 
termined from a simple Gaussian fit to the line, as a check 
on the results of the \ 2 minimization. 

3.2. C /f 3 P) XX1561 and 1657 multiplets 

The central portion of the 1561 A and 1657 A multiplets, 
arising from the 3 P ground term, were heavily saturated 
in the STIS data. However, this allowed us to examine the 
smaller variable red and blueshifted features in the spectra 
which are the signatures of infalling comets in the CS disk. 
In Figures 2 and 3, one can easily see absorption lying at 
higher rcdshift in the December 6 data that is not visible 
in the December 19 data. A model of the multiplet which 
contains only one velocity component is unable to repro- 
duce the absorption dip at 1561.5 A or the absorption at 
higher rcdshift in the December 6 data. 

Due to the fact that the multiplets were extremely sat- 
urated and that the multiple velocity components were 
blended together, we were not able to perform \ 2 min- 
imizations to find unique best values for the model pa- 
rameters. Fe II absorption features in our data set were 
analyzed and the velocities of their multiple velocity com- 
ponents found. Models of the C I A1561 multiplet were 
constructed, containing the same number of components 
as were found in the Fc II features, with roughly the same 
velocities. The parameters of the stable component at v 
= 20 km s _1 were set to those determined from the anal- 
ysis of the unsaturated 1613.376 A line. The remaining 
parameters were then adjusted to find the best model by 
eye. This model was compared to the 1657 A multiplct 
to confirm that the values found in this way were reason- 
able. A small constant value («2x 10~ 13 erg s _1 cm~ 2 
A -1 ) was subtracted from the 1657 A data to bring the 
baselines of the totally saturated features to zero. 

3.3. CI^D) W1931 and 1^63 lines 

These lines arise from the excited, metastable 1 D state 
of the ground configuration and have not previously been 
observed in spectra of (3 Pictoris. The 1 D level has a life- 
time of about 4000 s, lying about 10,000 cm -1 above the 
ground level; C atoms in this state may be produced dur- 
ing the photodissociation of CO in solar system comet co- 
mae (Tozzi, Feldman, & Festou 1998). The 1931 A lines 
were analyzed using the same model generating and x 2 
minimization procedures described in § 3.1; the atomic 
data used were from Hibbert et al. (1993). Although there 
may be some reason to suspect that the strong central ab- 
sorption features near 1931.05 A contain multiple veloc- 
ity components, since this behavior is seen in lines arising 
from excited levels of other atomic species, the x 2 mini- 
mizations indicated that stable unique solutions contain- 
ing more than one component in the central absorption 
features could not be found. Thus, best models were found 
using only one velocity component. However, the difficulty 
involved in modeling a saturated, blended line is such that 
there may well be other undetected components present, 
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so models with two components in the central absorption 
feature were compared to the data by eye. The 1463 A 
spectra, shown in Figure 5, were so noisy that they served 
only to roughly confirm the values found from the 1931 
A line. A small constant flux value 1 x 10~ 13 erg s _1 
cm~ 2 A -1 ) was subtracted from the 1463 A data to bring 
the baselines of the saturated lines to zero. 

3.4. CO Fourth Positive band system 

The (0-0), (1-0), and (2-0) bands of the Fourth Positive 
system of CO (A 1 H — X 1 E + ) did not vary between the 
two days of observation. The spectra were therefore aver- 
aged together to improve the S/N; the result appears in 
Figure 6. Only a single velocity component was observed. 
Note the detection of bands arising from 13 CO, the pertur- 
bation band e 3 £~ — X 1 H + (TO), not previously observed 
in spectra of (3 Pictoris, and the strong perturbation band 
d 3 A - 1^+ (5-0). 

Models were generated using wavelengths and oscillator 
strengths from Morton & Norcau (1994). Energies of the 
ground state levels were calculated using the Dunham co- 
efficients from Farrenq et al. (1991) and LTE assumed in 
order to determine the population of the rotational levels. 
The parameters of the model are the rotational excita- 
tion temperature, T, the Doppler broadening parameter, 
6, the column density of 12 CO in the ground vibrational 
state, 7V( 12 CO), the column density of 13 CO in the ground 
vibrational state, iV( 13 CO), and the velocity centroid, v. 
X 2 minimization was then performed on all three bands 
simultaneously. 

4. RESULTS 

4.1. C If 3 P) stable component 

The C I] A1613.376 line showed only one velocity com- 
ponent at 20 km s _1 , the systemic velocity of the star; the 
nominal uncertainty in the absolute wavelength calibra- 
tion of STIS leads to an error in velocity determinations 
of about 1 km s^ 1 (Sahu et al. 1999). This velocity and 
the fact that the line did not change between the two days 
of observation identifies the line as arising from stable gas. 
The results of the x 2 minimizations for all the features an- 
alyzed appear in Table 2. The column density for stable 
carbon in the ground level, 3 Po 5 determined from \ 2 min- 
imization was found to agree with that determined from 
the equivalent width of the line. Similarly the v and b 
determined from x 2 minimization agreed with the values 
found from fitting a simple Gaussian to the unsaturated 
line. 

Using a 3-cr upper limit on the column density of C I 
in the 3 P2 level (from the non-detection of the fine struc- 
ture line at 1614.5068 A) and assuming LTE, a firm upper 
limit of 100 K on the excitation temperature of 3 P carbon 
in the stable component was found. From analysis of the 
1561 A multiplet, discussed below, it was found that the 
excitation temperature of the stable component must be 
greater than about 50 K, or the multiplet could not be 
reasonably modeled. This range in temperature allows us 
to determine that the total column density of stable C I 
in the 3 P ground term is (2 — 4) x 10 16 cm~ 2 . 

This column density is more than an order of magni- 
tude larger than the total 3 P column density found by 
Jolly et al. (1998) from GHRS data taken in November 



1994, iV( 3 P) = 2 x 10 15 cur 2 . Either the abundance of 
stable C I has varied over the three years between obser- 
vations or, as is more likely, the difficult modeling of the 
heavily saturated C I A 1561 multiplet in the GHRS data 
led to an inaccurate column density. Models with column 
densities of 10 16 cm -2 in the stable component were com- 
pared to the 1994 GHRS spectrum of the 1561 A multiplet 
and found to fit the data equally well as models with the 
lower column density. We found the Doppler broadening 
parameter of the stable C I( 3 P) to be 1.3±0.5 km s" 1 . The 
previous work on the GHRS data found a Doppler broad- 
ening parameter of 4.2 km s _1 , which was much greater 
than the parameter found for the other atomic species and 
CO (Jolly et al. 1998). Our smaller b value is equal to that 
found for CO in our data; we thus do not see any excess 
kinetic energy in the motions of the C I atoms. 

4.2. C if^P) multiple velocity components 

The difficulty in modeling a saturated, blended multi- 
plet is formidable, and the values determined for the mul- 
tiple velocity components of C I( 3 P) from the 1561 A mul- 
tiplet arc not reliable. However, this analysis did confirm 
that models with column densities of 10 16 cm -2 in the sta- 
ble component could reasonably fit the 1561 A spectra. We 
can also conclude that the total column density of C I( 3 P) 
in the variable components is about 10 14 cm~ 2 on Decem- 
ber 6 and about 10 15 cm~ 2 on December 19, to within an 
order of magnitude or so. It is also clear that the variable 
C I features are generally better fit with higher excitation 
temperatures and larger b parameters, but these values are 
very poorly determined. 

4.3. Cl( l D) 

Using models with only one velocity component in the 
strong central absorption features, two components arising 
from the excited 1 D level of C I were found in the Decem- 
ber 6 data, only one in the December 19 data. The best 
models are overplotted on the spectra in Figures 4 and 5. 
Note that the central absorption features do not appear 
at the same velocity on the two days, that they also do 
not appear at the systemic velocity of the star, and that 
the column density of C I( X D) changes significantly be- 
tween the two days of observation. This would seem to 
indicate that there is no stable component for C I( 1 D). 
Models with two velocity components in the central ab- 
sorption feature, which were compared to the data by eye, 
indicate that there could be a "stable" component which 
has the same column density on both days. However, it 
must be at about 22-23 km s _1 , which is significantly dif- 
ferent from the systemic velocity of (3 Pictoris. These data 
are not able to conclusively determine the velocity struc- 
ture present in the C I( 1 D) gas; observation of unsaturated 
lines at similarly high resolution will be necessary. 

4.4. CO 

The CO bands showed no multiple velocity component 
structure and no change between the two days of ob- 
servation; the absorbing gas is thus entirely associated 
with the stable component at the systemic velocity of 
the star. The best model is shown overplotted on the 
data in Figure 6. The rotational excitation temperature 
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and the ratio of iV( 12 CO) to iV( 13 CO) show no signifi- 
cant change from the values found by Jolly et al. (1998). 
i?( 12 CO/ 13 CO) = 15±2 is quite small compared to typical 
values found in the ISM, e.g. R = 150 ± 27 for the diffuse 
clouds toward £ Ophiuchi (Shcffer et al. 1992). However, 
chemical fractionation at low kinetic temperatures can ex- 
plain this unusual ratio, as discussed below. The column 
densities of both 12 CO and 13 CO were found to be about 
1/3 the values found from the GHRS data. Note that the 
column density of 12 CO is about 50 times smaller than the 
column density of stable C I( 3 P). 

5. DISCUSSION 

Although CO in solar system comet comae is photodis- 
sociated by solar FUV photons, the CO in the (3 Pictoris 
CS disk is primarily destroyed by interstellar photons. The 
dissociation energy for CO is 11.1 cV; thus CO may only 
be dissociated by photons with wavelengths shortward of 
~ 1100 A. The type A5 star j3 Pictoris lacks the strong 
FUV emission lines created in the Sun's chromosphere; 
thus it emits very little flux in the FUV. The only source 
for CO-dissociating photons at (3 Pictoris is therefore the 
interstellar UV radiation field. Since CO photodissociates 
primarily through discrete line absorptions, self-shielding 
can have a strong effect on the abundance of CO in in- 
terstellar clouds (Van Dishoeck & Black 1988). However, 
this is unlikely to occur in the (3 Pictoris CS disk because 
of the very small transverse dimension of the disk. Thus, 
the photodissociation rate for 12 CO in the (3 Pictoris disk 
should be equal to the unshielded value and should also 
be equal to the value for 13 CO. Since 13 CO is not selec- 
tively dissociated in this situation, its high abundance was 
explained by the reaction 

i3 C++ i2 CO _ 13 CO+ 12 C + + 35K 

which favors the production of 13 CO at gas kinetic temper- 
atures below 35 K (Jolly et al. 1998). Assuming chemical 
equilibrium and that isotopic exchange is much more im- 
portant than photodissociation for both 12 CO and 13 CO 
(Shcffer et al. 1992), 



n( 12 CO) 
n( 13 CO) 



= exp 



-35 



12c 
13C 



15 ±2. 



The rotational excitation temperature of CO often does 
not accurately describe the gas kinetic temperature in 
diffuse environments (Wannier, Penprase, & Andersson 
1997). However, assuming a ( 12 C/ 13 C) ratio, we may use 
the above expression to estimate the true gas kinetic tem- 
perature of the CO. 

The average carbon isotopic ratio in the local ISM is 
~ 60 — 70 (Langer & Penzias 1993) and the typical solar 
system value found in comets is ( 12 C/ 13 C) = 89 (Je- 
witt et al. 1997). Using the range ( 12 C/ 13 C) = 89 - 60, 
we find that the gas kinetic temperature of the carbon 
monoxide is 20 K — 25 K, indicating that the CO gas is 
indeed colder than the stable C I gas. This suggests that 
the C I and CO are not located in the same regions of the 
disk. The assumption that isotopic exchange is more im- 
portant than photodissociation should be reasonable. For 
this assumption to apply, 



where T is the unshielded photodissociation rate, 2 x 
10~ 10 s- 1 (Van Dishoeck & Black 1988), k^ is the for- 
ward reaction rate, 6.8 x 10~ 10 cm 3 s _1 at 80 K (Smith & 
Adams 1980), and n( 12 C+) is the volume density of 12 C+. 
The average volume density of C, n( 12 C) = N( 3 P)/r c± 
20 cm~ 3 , assuming that the carbon extends over a distance 
r = 100 AU. Further assuming that n( 12 C+) > n( 12 C), 
as is likely, the right hand side of the above expression 
is greater than or equal to about 2 x 10~ 9 s _1 , an or- 
der of magnitude larger than T. However, a more com- 
plete treatment of the relationship between ( 12 CO/ 13 CO) 
and ( 12 C/ 13 C) takes into account photodissociation but 
requires exact knowledge of the density of 12 C + in the CS 
disk (Sheffer et al. 1992). 

The much larger column density of C I( 3 P) compared 
to that of CO leads us to believe that photodissociation of 
CO cannot be the only source of stable C I in the (3 Pictoris 
CS disk. Since there is no evidence that CO is produced by 
infalling comets (no variability or red and blucshiftcd fea- 
tures) , it has been postulated that the CO gas slowly evap- 
orates from the OEBs at several tens of AU from the star 
and is photodissociated to produce the stable C I (Lecave- 
lier des Etangs 1998). Obviously some portion of the C I 
gas must be produced directly from the FEBs (the portion 
giving rise to the variable red and blueshifted absorption); 
perhaps this C I gas is decelerated somehow and accumu- 
lates in the stable component before being destroyed by 
photoionization. In this scenario, the equilibrium column 
density of C I is 



n x N 



FEB 



r < k f 



exp 



-35 



n( 12 C+), 



N( 3 P) = 



where n is the mean number of infalling comets per year, 
~ 10 2 per year (Vidal-Madjar, Lecavelier des Etangs, & 
Fcrlet 1998), Nfeb is the total column density of C I gas 
produced by an infalling comet, and T is the photoion- 
ization rate for C I, 0.004 yr _1 . This expression, which 
assumes that C I atoms are lost only through photoion- 
ization by interstellar UV photons, gives N FEB ~ 10 11 
cm~ 2 , which should easily be produced by infalling comets 
like the ones giving rise to the variable C I( 3 P) compo- 
nents in our data. However, this is a rough treatment and 
the number of infalling comets per year varies significantly 
over time scales of a few years (Vidal-Madjar, Lecavelier 
des Etangs, & Ferlet 1998). Also, this treatment does not 
take into account reformation of C I by radiative recombi- 
nation; we cannot take this process into account properly 
without knowledge of the amount of C II in the CS disk 
and the electron density. 

Considering the C I( 1 D), when such atoms are produced 
by photodissociation of CO, O I( 1 D) atoms must be pro- 
duced also to conserve spin. The minimum total photon 
energy needed for this dissociation is 14.33 eV, correspond- 
ing to a threshold wavelength of 865 A. Since this thresh- 
old is below the Lyman limit, there are virtually no in- 
terstellar UV photons capable of producing C I( 1 D) by 
photodissociation of CO. Thus, the C I( 1 D) atoms in the 
(3 Pictoris CS disk cannot be produced by photodissocia- 
tion of CO by stellar or interstellar photons and must be 
produced by a collisional process involving ground state 
carbon atoms. Since the energy of the 1 D state relative to 
the ground state is high, the collisional process must be 
a very energetic one and therefore is likely to be closely 
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associated with the infalling bodies. Consequently, with a 
short C I( 1 D) lifetime, there would be no stable component 
in this gas; this behavior has not been previously seen in 
any constituent of the Pictoris CS disk. But this result 
is tentavive; examination of an unsaturated line at high 
resolution is needed to determine the velocity structure of 
the C I( 1 D) gas. 

6. CONCLUDING REMARKS 

The very high resolution, low scattered light contamina- 
tion, and good order separation of this STIS echelle data 
set has provided some clear advantages over previous ob- 
servations of Pictoris. The rotational lines of CO have 
been resolved, allowing for a much more precise determi- 
nation of the physical parameters of the gas. The column 
density of CO is N(CO) = (6.3 ± 0.3) x 10 14 cm" 2 and 
the ratio i?( 12 CO/ 13 CO) = 15 ± 2 is found. The absence 
of transient red or blueshifted components in the high res- 
olution CO spectra supports the suggestion that this gas 
evaporates from cometary bodies orbiting far (several tens 
of AU) from the star. But the fact that the column den- 
sity of CO is only about 2% of the total column density 
of C I in the 3 P ground term implies that photodissocia- 
tion of this CO is not the primary source for C I gas. It 
could perhaps be produced directly from infalling comets 
close to the star, but the mechanism by which it comes to 
zero velocity relative to the star and accumulates before 
being photoionized is unclear (although see Lagrange et 
al. (1998)). The C I( 1 D) gas may not have a stable com- 
ponent at 20 km s ; this unique species could prove to 
be a valuable tracer of FEB activity in the Pictoris CS 
disk. 

Despite the advantages of this data set, our lack of suc- 



cess in modeling the heavily saturated C I multiplcts indi- 
cates that in order to really determine the characteristics 
of the variable components of the C I gas, we need to ob- 
serve an unsaturated line or multiplet, with an oscillator 
strength between that of the 1561 A multiplet and that of 
the spin-forbidden 1613.376 A line. A number of suitable 
multiplcts and lines lie in the FUV, shortward of ~ 1300 A. 
Also, an unsaturated line arising from the 1 D level would 
allow us to confirm the velocity structure in this gas and 
to determine if the velocities of the X D gas components 
correspond with any of the velocities of the variable com- 
ponents in the 3 P gas. Three likely lines lie between 1311 
A and 1359 A. Measurement of the densities of C II and 
O I would greatly help to unravel the carbon chemistry of 
the Pictoris disk. Again, potentially useful multiplcts of 
these species lie in the FUV below 1340 A. Thus, although 
this data set has vastly increased our knowledge about the 
important species C I and CO in the Pictoris disk, our 
understanding would probably benefit greatly from inves- 
tigation of Pictoris at shorter ultraviolet wavelengths. 
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is operated by the Association of Universities for Research 
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grant GO-07512.01-96A from the Space Telescope Science 
Institute. 



REFERENCES 



Beust, H., Lagrangc-Hcnri, A. M., Vidal-Madjar, A., & Ferlet, R. 

1990, A&A, 236, 202 
Draine, B. T. 1978, ApJS, 36, 595 

Farrenq, R., Guelachvili, C, Sauval, A. J., Grevesse, N., & Farmer, 

C. B. 1991, J. Molec. Spectrosc, 149, 375 
Hibbert, A., Biemont, E., Godcfroid, M., & Vaeck, N. 1993, A&AS, 

99 179 

Jewitt, D. C, Matthews, H. E., Owen, T., & Meier, R. 1997, 

Science, 278, 90 
Jolly, A. et al. 1998, A&A, 329, 1028 
Lagrange, A. M. et al. 1998, A&A, 330, 1091 
Langer, W. D., & Penzias, A. A. 1993, ApJ, 408, 539 
Lecavclicr des Etangs, A., Vidal-Madjar, A., & Ferlet, R. 1996, 

A&A, 307, 542 
Lecavelier des Etangs, A. 1998, A&A, 337, 501 



Morton, D. C. 1991, ApJS, 77, 119 
Morton, D. C. & Noreau, L. 1994, ApJS, 95, 301 
Sahu, K. et al. 1999, STIS Instrument Handbook v3.0 (Baltimore: 
STScI) 

Sheffer, Y., Federman, S. R., Lambert, D. L., & Cardclli, J. A. 

1992, ApJ, 397, 482 
Smith, B. A., & Terrile, R. J. 1984, Science, 226, 1421 
Smith, D., & Adams, N. G. 1980, ApJ, 242, 424 
Tozzi, G. P., Feldman, P. D., & Festou, M. C. 1998, A&A, 330, 753 
Van Dishoeck, E. F. & Black, J. H. 1988, ApJ, 334, 771 
Vidal-Madjar, A. et al. 1994, A&A, 290, 245 

Vidal-Madjar, A., Lecavelier des Etangs, A. L., & Ferlet, R. 1998, 

Planet. Space Sci., 46, 629 
Wannier, P., Penprase, B. E., & Andersson, B-G 1997, ApJ, 487, 

L165 



6 



Roberge et al. 



Table 1 



HST-STIS OBSERVATION LOG 



Observation 


Grating 


Slit 


Spectral Range 


Start Time 


Exposure Time 


ID# 




(arcsec) 


(A) 


(UT) 


(s) 


1997 Dec 6 














o4g001010 


E140H 


0.2 x 0.09 


1461 - 


1663 


07:17:58 


900.0 


o4g001020 


E230H 


0.1 x 0.03 


1874 - 


2146 


07:40:26 


80.0 


o4g001030 


E230H 


0.1 x 0.09 


1624- 


1896 


07:48:20 


678.8 


o4g001040 


E230H 


31 x 0.05 


2624- 


2895 


08:48:28 


360.0 


o4g001050 


E230H 


31 x 0.05 


2374- 


2646 


09:02:20 


360.0 


o4g001060 


E230H 


31 x 0.05 


2124 - 


2396 


09:16:12 


288.0 


1997 Dec 19 














o4g002010 


E140H 


0.2 x 0.09 


1461 - 


1663 


19:51:29 


900.0 


o4g002020 


E230H 


0.1 x 0.03 


1874- 


2146 


20:13:57 


80.0 


o4g002030 


E230H 


0.1 x 0.09 


1624- 


1896 


20:21:51 


678.8 


o4g002040 


E230H 


31 x 0.05 


2624- 


2895 


21:25:45 


360.0 


o4g002050 


E230H 


31 x 0.05 


2374- 


2646 


21:39:37 


360.0 


o4g002060 


E230H 


31 x 0.05 


2124- 


2396 


21:53:29 


288.0 
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Table 2 
Results of Modeling 





Velocity 


Column density 


Doppler broadening parameter 


Temperature 




v (km s ) 


AT ( mn — 2 \ 


o ^Km s ) 


1 (KJ 


Stable C I( 3 P ) 












20 ± 1 


(6±1) x 10 15 


1.3 ±0.5 




C I( 3 P) 










December 6 


20 ± 1 


(2-4) x 10 16 


1.3 ±0.5 


50 - 100 




26 


10 14 


4 


10 2 




41 


10 14 


10 


10 3 




57 


10 13 


12 


10 3 


December 19 


11 


10 13 


5 


10 1 




20 ± 1 


(2-4) x 10 16 


1.3 ±0.5 


50 - 100 




27 


10 15 


5 


10 2 


CI^D) 










December 6 a 


zo ± 1 


^4 ± 6) X 1U 


(A 4- 1 "\ 

(4± 1) 






49 ±2 


(4 ±2) x 10 12 


(4 ±2) 




December 6 b 


22 


4 x 10 13 


2 






28 


8 x 10 12 


3 






48 


6 x 10 12 


3 




December 19 a 


26 ± 1 


(1.2 ±0.7) x 10 14 


4.4 ±0.8 




December 19 b 


23 


4 x 10 13 


2 






29 


3 x 10 13 


3 




12 CO 












20 ± 1 


(6.3 ±0.3) x 10 14 


1.3 ±0.1 


15.8 ±0.6 


13 co 












20 ± 1 


(4.3 ±0.4) x 10 13 


1.3 ±0.1 


15.8 ±0.6 



a Results from x 2 minimizations using model with one velocity component in the central absorption 
features 

b Results from chi-by-eye using model with two velocity components in the central absorption 
features. 
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Fig. 1. — The C I] A1613.376 line, with both days of data averaged together. Positions of the two other non-detected fine structure lines in 
the multiplet are indicated as well. The spectrum has not been rebinncd or smoothed. 
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2.0 r 




1560.0 1560.5 1561.0 1561.5 1562.0 1562.5 

Wavelength (A) 



Fig. 2. — The C I A1561 multiplet, with both days of data shown. The data have not been rebinned or smoothed. The continua around the 
features were fit with sixth-degree polynomials and the best fitting models based on the C I( 3 P) parameters in Table 2 are shown overplotted 
on the data. 
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t i i i i i i i i i i y__, i i i i i i i i i i i i i i d 

1656.0 1656.5 1657.0 1657.5 1658.0 1658.5 

Wavelength (A) 



Fig. 3. — The C I A1657 multiplet, with both days of data shown. The data have not been rebinned or smoothed. The continua around the 
features were fit with sixth-degree polynomials and the best fitting models based on the C I( 3 P) parameters in Table 2 are shown overplotted 
on the data. 



CO and C I in the (3 Pictoris CS disk 



11 



Velocity (km s ') 

-40 -20 20 40 60 80 100 120 

2.0 r i ' i ' i ' i ' i ' i ' i ' i ' r~ 




1930.6 1930.8 1931.0 1931.2 1931.4 1931.6 

Wavelength (A) 

Fig. 4. — The C I A1931 line, with both days of data shown. The data have not been rebinned or smoothed. Heliocentric velocity appears 
on the upper x-axis. The continua around the features were fit with fourth-degree polynomials and the best fitting models based on the 
C I( X D) parameters in Table 2 are shown overplotted on the data. 



12 



Roberge et al. 



Velocity (km s ') 

-80 -60 -40 -20 20 40 60 80 100 120 14C 




1463.0 1463.2 1463.4 1463.6 1463.8 1464.0 

Wavelength (A) 



Fig. 5. — The C I A1463 line, with both days of data shown. The data have not been rebinned or smoothed. Heliocentric velocity appears 
on the upper x-axis. The continua around the features were fit with fourth-degree polynomials and the best fitting models based on the 
C I( X D) parameters in Table 2 are shown overplotted on the data. 
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1510.5 
Wavelength (A) 



CO A-X (2-0) 




1478.0 1478.5 
Wavelength (A) 



Fig. 6. — The CO A — X (0-0), (1-0), and (2-0) bands, with both days of data averaged together. The data have not been rebinned or 
smoothed. The total model based on the CO parameters in Table 2 is overplottcd as a dashed line; the contribution from 13 CO alone is 
ovcrplottcd as a solid line. 



